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We show how statistical thermodynamics can be formulated in situations in which thermodynamics 
applies, while equilibrium statistical mechanics does not. A typical case is, in the words of Landau 
and Lifshitz, that of partial (or incomplete) equilibrium. One has a system of interest in equilibrium 
with the environment, and measures one of its quantities, for example its specific heat, by raising the 
temperature of the environment. However, within the observation time the global system settles down 
to a state of apparent equilibrium, so that the measured value of the specific heat is different from the 
equilibrium one. In such cases formula for quantities such as the effective specific heat exist, which are 
provided by Fluctuation Dissipation theory. However, what is lacking is a proof that internal energy 
exists, i.e., that the fundamental differential form 5Q — 5W (difference between the heat absorbed by the 
system and the work performed by it) is closed. Here we show how the coefficients of the fundamental 
form can be expressed in such a way that the closure property of the form becomes manifest, so 
that the first principle is proven. We then show that the second principle too follows, and indeed 
as a consequence of microscopic time-reversibility. The treatment is given in a classical Hamiltonian 
setting. One has a global time- independent Hamiltonian system constituted by the system of interest 
and two auxiliary ones controlling temperature and pressure, and the occurring of a process due to a 
change in the thermodynamic parameters is implemented by a suitable choice of the measure for the 
initial data. 



I. INTRODUCTION 

There are situations of metaequilibrium or metastabil- 
ity, in which the laws of phenomenological thermodynam- 
ics apply, whereas equilibrium statistical thermodynam- 
ics does not. Here we show how at least a first step can be 
made in providing a formulation of statistical thermody- 
namics for such cases, working in a classical Hamiltonian 
setting. 

A paradigmatic case is that of partial (or incomplete) 
equilibrium, as Landau and Lifshitz call it (see 1, sec. 4). 
Macroscopically one has a system characterized by a phe- 
nomenological internal energy, say U{T,p) if its thermo- 
dynamic state is defined by temperature T and pressure 
p. Microscopically one has a corresponding dynamical 
system with a set of degrees of freedom, a Hamiltonian 
H and a given family of invariant probability measures, 
for example the generalized Gibbs measure with parame- 
ters T and p. At equilibrium one makes the identification 



U 



<^)t,p 



where {■) denotes mean with respect to the invariant 
generalized Gibbs measure for the given values of T and 
p. In the case of partial (or incomplete) equilibrium, how- 
ever, up to a given observation time r only a fraction of 
the degrees of freedom is actually able to react to changes 
of the thermodynamic parameters. 
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Thus, considering for example the case of the specific 
heat at constant volume Cy , this is a well defined thermo- 
dynamic quantity which has a concretely measured value, 
but this value does not agree with the one predicted by 
equilibrium statistical mechanics. A typical case, men- 
tioned by Landau and Lifshitz, is "the partial equilib- 
rium of a mixture of several substances which interact 
chemically. Owing to the relative slowness of chemical 
processes, the equilibrium connected with the motion of 
molecules occurs, in general, considerably sooner than 
the equilibrium connected with the interchange of atoms 
between the molecules, i.e., connected with the compo- 
sition of the mixture". For example, in a system com- 
posed of a mole of graphite C and a mole of O2, at 
room temperature combustion does not occur, and the 
measured heat capacity at constant volume is simply 
the sum of the specific heats of the two components, 
i.e., 1R+ (5/2)i? = {7/2)R, whereas at equilibrium the 
heat capacity is that of oire mole of carbon dioxide CO2 , 
namely, (5/2)i?. Another situation presenting some anal- 
ogy with partial equilibrium is that of the "time depen- 
dent specific heat" of Birge and Nagel [2]. 

On the other hand, in all the cases mentioned there 
exists a theoretical formula for Cy which explicitly de- 
pends on the observation time, namely, the formula pro- 
vided by Fluctuatioir Dissipation (FD) theory, which is 
thus a good candidate for a theoretical counterpart to the 
phenomenological value. More in general, in situations of 
metaequilibrium one has available the differential forms 
SQ and 6W of the heat absorbed by the system and of the 
work performed by it, with coefficieirts which have well 
defined phenomenological values. Analogously, the corre- 
sponding theoretical coefficients provided by FD theory 
have well defined values too, if a suitable dynamical re- 



laxation did occur, albeit possibly a partial one. 

However, there remains open the problem of proving 
that a consistent statistical thermodynamics can be for- 
mulated. At a phenomenological level, the internal en- 
ergy can be defined by integration of the available differ- 
ential form 6Q — SW, which is closed in virtue of the first 
principle. At a theoretical level, instead, such a property 
should be proved, but a proof is apparently lacking. 

In this paper we show that, in situations of metae- 
quilibrium, the formulae for the coefficients of the fun- 
damental form SQ — SW provided by FD theory can be 
expressed in such a way that the closure property of the 
form becomes manifest. This proves the first principle. 
Moreover, the second principle too is proven, as a con- 
sequence of microscopic time-reversibility. So one would 
have a complete statistical thermodynamics if one could 
also prove the validity of the zeroth principle (transitivity 
of metaequilibrium) , or of a suitable generalization of it. 
In any case one has available a candidate for the internal 
energy, namely, the state function U which is obtained by 
integrating the form SQ — SW, with coefficients provided 
by FD theory. However, and this is perhaps the more 
relevant contribution of this paper, it turns out that in 
general such an internal energy U cannot be expressed as 
the mean of the Hamiltonian with respect to the invariant 
measure considered. 

The main point thus consists in introducing a suitable 
definition for the statistical mechanical analogues of Q 
and W, as quantities parametrically dependent on the 
observation time r (which we think of as fixed once for 
all, and will be mostly left understood). To this end 
we essentially follow the procedure used in FD theory, 
which in turn mimics the procedure of phenomenologi- 
cal thermodynamics, inasmuch as it takes into account 
only the fraction of energy which is actually exchanged 
with the measurement apparatus, when an increment is 
given to the thermodynamic parameters T and p of the 
latter. It is left to the dynamics itself, which involves the 
interactions between the system of interest and the mea- 
surement apparatus, to determine the amount of energy 
which is actually exchanged. It will be seen, however, 
that the definitions given here are not exactly those of 
FD theory. Some criticism was already raised against 
the latter [3J. The present definitions are in our opinion 
somehow simpler, and more suited to the aims of ther- 
modynamics. 

In order to implement the idea of mimicking the proce- 
dure used in phenomenological thermodynamics, we have 
first of all to consider the system of interest as coupled 
with two auxiliary ones, a thermostat which controls the 
temperature and a barostat which controls the pressure. 
The global system is studied as an isolated, time inde- 
pendent, Hamiltonian system. The existence of differ- 
ent states corresponding to different values of the macro- 
scopic thermodynamic parameters T and p is taken into 
account by introducing different choices for the probabil- 
ity measure of the initial data in the global phase space. 
This way of handling the problem will make the proofs of 



the first and of the second principles particularly simple. 

In section |ll] the analogue of the fundamental form is 
defined and the first principle is proven. The second prin- 
ciple is proven in section iHll To this end, making use of 
microscopic time-reversibility, we preliminarily establish 
formulas of FD theory type, which express the coefficients 
of the form SQ in terms of equilibrium correlations of 
suitable functions. From this, in particular, the property 
Cp > follows, which amounts to the classical formula- 
tion of the second principle in the form of Clausius. In 
section HV] it is shown how the thermodynamic tempera- 
ture of the system of interest coincides with that of the 
thermostat if a suitable dynamical relaxation takes place. 
Further discussions are given in the conclusive section. 

In an appendix it is discussed how the present results 
may be of interest for the Fermi-Pasta-Ulam problem, 
the analogy of which with glassy-like systems was first 
proposed by a group of people about Parisi in the year 
1982 ^] (see also[^ and the review i6|). 



II. THE MODEL, MICROSCOPIC ANALOGUE 

OF THE FUNDAMENTAL FORM, AND 

DEDUCTION OF THE FIRST PRINCIPLE 

As usual (see 0, or |^, page 131), we take the system 
of interest (say, a fluid) enclosed in a cylindrical vessel 
one basis of which is fixed, while the other one is a mov- 
able piston. The distance of the piston from the fixed 
wall will be denoted by /, so that to any distance / there 
corresponds a volume T^ = ^/ of the system of interest, 
A being the area of the piston. Through the two bases 
of the cylinder the system of interest can exchange en- 
ergy with two other systems, a thermostat which acts as 
a source of heat (through the fixed basis), and a barostat 
which acts as a source of purely mechanical, or adiabatic, 
work (through the movable piston). It is well known [9| 
that in dynamical terms the adiabaticity condition is very 
well satisfied up to extremely long times if the piston is 
modeled as a single rigid particle of macroscopic mass. 
The two auxiliary systems are thought of as having no 
direct mutual interaction. 

The global system (system of interest plus sources of 
heat and of work) is dealt with in a classical Hamiltonian 
setting, with a given time independent Hamiltonian and 
a fixed global geometry, the position of the piston being 
considered as one of the configurational coordinates of 
the global system. One thus has a global phase space 
with canonical coordinates 

z = {x,x',x",l,pi) 

where x are the canonical coordinates of the system of 
interest, x' those of the heat source, x" those of the work 
source, I the distance of the piston from the fixed wall, 
and pi the corresponding conjugate momentum. The to- 
tal Hamiltonian is decomposed as 

Htot{z) = H{X, I) + Hti^cAx') + i?bar(a;", I) + Hintcr{z) ■ 



Here H is the partial Hamiltonian of the system of inter- 
est, while -father 1 and i?bar are the partial Hamiltonians of 
the sources, the labels ther and bar standing for thermo- 
stat and barostat respectively. Finally there appears an 
interaction Hamiltonian, which should contain a term of 
the formp^/2Af for the kinetic energy of the piston, dealt 
with as a single rigid particle of mass M of macroscopic 
size. As usual it is understood that, for large systems, 
the interaction Hamiltonian, although playing an essen- 
tial dynamical role in allowing for energy exchanges, can 
be neglected in estimating mean values of the quantities 
of interest. Indeed such quantities are proportional to 
the volume of the system of interest, whereas the inter- 
action terms are proportional to the area of the basis 
of the cylinder; furthermore, motions are considered in 
which the piston moves very slowly. The flow in the 
global phase space, induced by the global Hamiltonian 
iJtot, will be denoted by $*. 

Our aim is to estimate the energy exchanges between 
the system of interest and each of the sources (namely, 
the heat Q absorbed by the system and the work W per- 
formed by it). This involves first of all the dynamics 
of the global system, which induces corresponding ex- 
changes for any initial datum z, and then the specifica- 
tion of the measure over the initial data which should be 
assigned in order to estimate mean values. 

For what concerns the dynamics, along any solution of 
the equations of motion the total energy is conserved, 
so that, neglecting the contribution of the interaction 
Hamiltonian. one has energy conservation in the form 

[H + i/thcr + i^bar] ($"^) ^ [H + i/thcr + i^bar] (z) , 



or equivalently 



[Ai/ + AiJthcr + AHbar] W - , 



where we have introduced the notation 



AF{z) = A-'Fiz) =^ F($^2) - F{z) 



(1) 



(2) 



for the change of any dynamical variable F induced, at 
the observation time r, by the time evolution <i>* of the 
global system, as a function of the initial point z of phase 
space. Relation ([Ij will finally lead to a deduction of the 
first principle for the energy exchanges. This, however, 
requires to previously define the measure on the global 
phase space that should be used in estimating mean val- 
ues, or even typicality of an orbit. 

Clearly, invariant measures correspond to equilibrium 
states. Indeed the mean values of dynamical variables 
not explicitly depending on time, computed with respect 
to an invariant measure, turn out to be time indepen- 
dent. Thus the mean values of the three energies H, 
i?ther and i?bar do uot chaugc with time, i.e., there are 
no energy exchanges, notwithstanding the fact that no 
one of the three partial Hamiltonians is an integral of 
motion. Choosing a measure which is not invariant, leads 
instead to a thermodynamic process. We are interested 



in processes induced by an increment of the thermody- 
namic parameters of the sources, with respect to those 
characterizing an invariant measure. 

To this end we consider a family of measures relative to 
the global system. The family depends on two "external" 
thermodynamic parameters relative to the sources (say 
temperature T for the heat source and pressure p for the 
work source), and furthermore on at least two "internal" 
parameters, say a = (011,02), relative to the system of 
interest. So the family will have densities of the form 

p{T,p,a;z) . 

We then assume that, for any pair T,p of external 
thermodynamic parameters, there exists an equilibrium 
state of the global system (that might be called refer- 
ence state), which corresponds to some definite values 
a — a(T, p) of the internal parameters a for which the 
measure is invariant. In particular the invariant mea- 
sure could be a global Gibbs one. This is the typical 
case one considers in discussing partial equilibrium. The 
more general case considered above is of interest when 
discussing metastability. 

Given one such equilibrium state with density 



p{T,p,a;z) 



(3) 



we then introduce the state in which the external pa- 
rameters have been incremented, while the values of the 
internal parameters are kept constant. In general such 
a measure will not be invariant, i.e., will correspond to 
a nonequilibrium situation, leading to a process. The 
corresponding measure will have a density 



p{T + dT, p + dp,a;z) 



(4) 



Notice that, in the present model, a change of the exter- 
nal parameters can be obtained without any alteration 
of the global geometry of the system, because a change 
in the pressure p of the barostat can be produced by a 
suitable change of its temperature, which in turn can be 
described through a change of the measure on the initial 
data. 

Mean values with respect to the noninvariant measure 
with density (j4]) will be denoted by (•) , while those 
referring to the invariant measure with density ([3]) might 
be denoted by ( • ) , and instead will be simply denoted 



(F)"'^ dof r p^^^ ^^p + dT,p + dp, a; z) dz , 



(5) 



(F)'i^' / F{z)piT,p,a;z)dz 



Let us now take as a measure for the initial data of 
the global system the noninvariant one with density ((4]). 
For any initial datum z, which evolves to <I>^z, there 
remain defined an increment of the energy of the ther- 
mostat and an increment of the energy of the barostat. 



and this naturally leads to microscopic definitions of the 
heat Q absorbed by the system of interest and of the 
work W performed by it. Indeed, by analogy with the 
standard equilibrium procedure, for any given dynamical 
variable F[z) one considers the increment corresponding 
to each initial datum, and then takes the mean value over 
the data distributed according to the given (nonequilib- 
rium) measure. As usual, one might presume that such 
mean values are good representatives of the values actu- 
ally taken in correspondence to single typical orbits. 

So, recalling the meaning ^ of A_F, we introduce the 
definitions which constitute the key point of the whole 
paper, namely, 

- I [Hthcr($"z) - i/thor(^)] p{T + dT,p + dp, a; z) dz , 

(6) 

We can now come to the first principle. From the def- 
initions ([nj of Q and W , and from energy conservation 
dl]), one gets 



T4^ = (Ai/(z)) 



ncq 



I AH{z) p{T + dT, p + dp, S; z) Az 



(7) 



Then, by expanding p to first order in AT and dp, dif- 
ferential forms 5Q and 5W remain thus defined. So the 
analogue of the fundamental form 5Q — SW turns out to 
be expressed in terms of the Hamiltonian H of the sys- 
tem of interest only (although still in terms of a measure 
involving the two sources), being given by 



6W^ci{T,p)dT + C2{T,p)dp 



with coefficients 



ci = / AH{z) —p{T,p,a;z)dz , 
C2 = / AH{z) —p{T,p,a]z)dz . 



(8) 



(9) 



Proving the first principle amounts to proving that the 
differential form SQ — SW is closed. On the other hand, 
the variables T and p enter the integrals at the right hand 
side of © only through the factor -Mpp or the factor -S-p. 
Thus, assuming that the density p is smooth enough in 
the parameters, one has 



dci 
dp 



dc2 
'df 



which is the analogue of a Maxwell relation. So the form 
6Q — SW is closed, and there exists (at least locally) a 
function U = U{T,p) such that SQ - SW = dU. So the 
internal energy can be defined by integrating the differ- 
ential form ([5]). The first principle is thus proven. Notice 



that the proof is very simple due to the fact that we 
are considering a closed model, in which the changes in 
the external parameters involve only changes in the mea- 
sures, and not in the dynamics (at variance with what 
occurs in standard FD theory). 

We close this section with two remarks. The first one 
concerns the fact that the internal energy U cannot be 
expressed as the mean of the energy H with respect to 
the invariant measure. For example, the state function 
U defined by 

U= H{z) p{T,p,a;z)dz 



(10) 



produces a differential form 



dil = ci (T, p) dT + C2 (T, p) dp , 



with coefficients 

_ dcf 
Cl = 



ii{z) —p{T,p,a;z)dz 



C2 = / H{z) —p{T,p,a;z)dz 



(11) 



which, at variance with ([9]) , involve integrals of H rather 
than of AH. 

The second remark concerns the relation between the 
dynamical increment of internal energy 



A^C/ =^ Q 



W 



(12) 



defined through ([7]) ^ in correspondence to increments 
dT, dp of the parameters - and the increment of the equi- 
librium internal energy which is obtained by increment- 
ing the parameters T and p in the equilibrium measure. 
This relation is obtained by remarking that the definition 
([7]); which involves the time evolved H{^'^z) of the dy- 
namical variable H while keeping fixed the measure, can 
also be dually expressed by letting the measure evolve 
with time, while keeping fixed the dynamical variable H: 



A^U= I H{z) p{T + dT,p + dp,a;^-^z) 



- p[T + dT, p + dp, a; z) 



(13) 



Thus, if the density p tends (in weak sense) to the fi- 
nal equilibrium density, with the internal parameters a 
adapted to the new external ones, and if in addition it fac- 
tors into an external part and an internal one (as should 
be assumed of the initial density too), then one should 



have A^U ^ (h) 



T+dT,p+dp 



(H) 



T,p'' 



I.e. 



{H)^^+A^U^{H) 



T+dT,p+dp ' 



where ( • )^ denotes mean value over the phase space 
of the system of interest, with respect to the generalized 
Gibbs measure at values p, T of the parameters. 



So the time-dependent family /S^U considered here 
provides an interpolation (-ff)^ + A'^C/ between the 

initial equilibrium value iH\„ and the final one 



T+dT,p+dp 



of the internal energy. However, it is left 
to the dynamics to decide whether, and at which time, 
will the final equilibrium value be actually attained. 



III. THE SECOND PRINCIPLE AS A 

CONSEQUENCE OF MICROSCOPIC 

TIME REVERSIBILITY. 

The result proved in the previous section, according to 
which in general the internal energy of metaequilibrium 
thermodynamics cannot be expressed as the mean of the 
Hamiltonian with respect to the invariant measure, is 
the one to which we attach a particular significance, as 
will be discussed in the appendix. In this section and 
in the next one, however, we continue the discussion of 
metaequilibrium thermodynamics. Here we show how 
the second principle too is obtained. To this end we 
will show that the coefficients of the fundamental form 
6Q — SW can be expressed in terms of suitable correlation 
functions, as in FD theory. This fact is obtained if the 
further assumptions are made that i) the dynamics is 
time reversible, and thatii) the reservoirs are described 
by Gibbs measures. The latter condition can actually be 
weakened, as will be mentioned later. 

So, let us recall the time-reversibility property. This 
property amounts to requiring that there exists a map- 
ping V of the phase space onto itself (inverting the sign 
of the velocities of all particles that constitute the global 
system) with the properties 



V" =1 



(7?$^)2=x, 



where T is the identity map. Obviously the invariance 
will be required to hold also with respect to the con- 
sidered measure. It is well known that the reversibility 
property is satisfied by essentially all systems of interest, 
namely, those with Hamiltonians even in the momenta of 
the particles. So in particular in our case we will have 



HtotiVz) 

HtheriVz) 



Htotiz) , H{Vz) 
- Hthci{z) , 



H{z), 



Concerning the family of measures for the initial data, 
we will assume that the invariant measures have the form 






-/3bar-H"b 



r,syst 



(a;z) 



(14) 



where p^'^^*^ refers to the system of interest, c(/3,/3bar) is 
the familiar normalization factor involving the partition 
functions of the reservoirs, while /? and /3bar are Lagrange 
multipliers which determine the temperature T of the 
thermostat and the temperature of the barostat (and so 
implicitly its pressure too). 



Our aim is now to find a suitable expression for W 
and thus for Q. First of all the equilibrium pressure p is 
defined as usual (see Q) as the equilibrium mean value of 

the corresponding dynamical variable P{z) = f^(^) = 
-j^iz). So the pressure p is given by 

Thus the work performed by the system of interest on 
the barostat should be given by 



W^p{AV) 



ncq 



(15) 



So the heat Q absorbed by the system of interest, Q — 
(Ai/) -I- W (see ([7])), will eventually be expressed in 
terms of dynamical variables of the system of interest, as 

Q=(A(i/+p^)rv 

We thus turn to investigating the expression for the 
specific heat. We even perform the more general compu- 
tation which gives the differential form SQ in terms of d/3 
and d/3bar- Such a form is obtained from Q through the 
same procedure used in the previous section for the fun- 
damental form SQ — SW, i.e., by expanding the nonequi- 
librium density to first order in d/3 and d/3bar- So one 
gets 

SQ = gi(/3,/3bar)d/3-|-92(/3,/3bar)d/3bar , (16) 

where the coefficients are given by 



A{H+pV)iz) — p(/?,/3bar,a;^)dz , 
A{H+pV)iz) p(/?,/3ba..,a;z)dz 

O/Sbar 



(17) 



One now has the following 
Lemma: Let the global Hamiltonian system be time- 
reversible, and let the invariant measure have a density 
of the form p^P- Then, for any dynamical variable F 
having the property 



F{Pz) - F(z) , 



J AF ^pAz = ~\{^F A{H + pV)) 



(18) 



(19) 



The proof, for example for the first formula, is as fol- 
lows. Due to the exponential form of p with respect to 
the reservoirs one has 

/ '= ^(AF) = - j AFiz)Ht^,,{z)p{z)dz . 

Notice that, in perfornnng the partial differentiation, one 
further term would appear, which contains the derivative 



of the normalization factor c(/3,/3bar)- However, being 
proportional to (AFY such a further term vanishes, in 
virtue of the invariance of the considered measure. In- 
troduce then the change of variables 

z = T'^^y 

which is canonical, so that dz = dy. Use (fT5|. which 
implies 

AF{z) = - AF{y) , 

and use also the invariance of the measure, p(z) = p{y)- 
Thus, calling again by z the dummy variable y one gets 

/ = J AFiz)Hti,eriV<P^z)p{z)dz . 

Finally, the result follows by taking for / the semisum of 
the original expression and of the last one, using 



temperature of the thermostat will be provided in the 
next section. 

Finally, the fact that the thermodynamic quantities 
may attain the corresponding equilibrium values if a suit- 
able relaxation takes place, is exhibited by expressing the 
thermodynamic coefficients in terms of time correlation 
functions. For example, in the case of the specific heat 
Cp given by (j2ip . by expanding the square 

[A{H+pV)iz)f = [{H+pV){^^z) - {H+pV){z)f 



and adding and subtracting the term (^H 
immediately gets 



pvy 



G„ — 



ksT^ 



^H+pV - (m+pvir) 



(22) 



(F(<i>'^z)F(z)) — (F) its time-autocorrelation function. 



where ap denotes the variance of F and 



AHi 



the 



^A{H+pV) . 



The second formula is proved analogously, using 
AHh^,=pAV. 

Thus the coefficients gi, 52 of the form dQ turn out 
to be expressed in terms of correlation functions with 
respect to the invariant measure: 



qi = -l{[AiH+pVf) 
q2 = ^{A{H+pV)AV) . 



(20) 



In particular, using ^ = —il/ikBT'^)^ one finds that 



the specific heat Cp at constant /3ba 
is manifestly positive. 



(i.e., at constant p). 



Gt, 



1 



2kBT- 



-{[A{H+pV)f) >0 



(21) 



Thus the second principle in the Clausius form (heat 
spontaneously flows from hot to cold bodies) is proven. 
By the way, one is naturally led to guess that the second 
principle may not hold if microscopic time-reversibility 
were not satisfied. 

The positiveness of the specific heat was proven above, 
making use of the assumption that the marginal measure 
for the environment be of Gibbs type. However, if one 
looks at the proof one realizes that the result also holds if 
any assumption is made which guarantees that the partial 
derivatives of the density p have negative definite sign. 
For example, a sufficient condition is that, in formula 
(fT4| for the density, instead of the exponential exp{~(3H) 
there appears a factor f{/3H), with any function / having 
the property /' < 0. 

Having thus proven positiveness of the specific heat, 
even in such a more general frame, the existence of an in- 
tegrating denominator for the form SQ then follows from 
arguments of Caratheodory type. A connection with the 



IV. THE TEMPERATURE OF THE 

THERMOSTAT AS AN INTEGRATING 

DENOMINATOR FOR SQ 

Thus, having proven the second principle, we know 
that an integrating denominator for the differential form 
6Q exists. One may ask however whether such an inte- 
grating denominator coincides with the temperature 1//3 
of the thermostat. We will show now that this is actually 
the case, at least if the observation time r is sufficiently 
large to guarantee that a suitable dynamical decorrela- 
tion did occur. 

So, let us consider the form 



pSQ = /?(7i d/3 + I3q2 d/?ba 



(23) 



where the coefficients Qi are defined by ([T7|) . and let us 
investigate whether it is closed. This requires to calculate 

Qg dT • ^^ ^^ ^^^ '"^^^ '^^ ^^^ internal energy, 

again one finds that the two terms involving the second 
derivatives of p cancel, so that one remains with 



d{(3qi) 9(/3<72 



9/3bar d/3 

dp 



P- 



j AV{z] 



- j A{H+pV){z) 

So, using the Lemma one gets 
-d{Pqi) d{Pq2)- 



d_ 

8(3 
d 



a/3, 



pdz 
pdz 



(24) 



bar 



9/3b 



dp 



= --{AV AiH+pV)) 



dp 



bar 



Thus the form is closed if one has 
dp p 

9/3bar " " ^ 



(25) 



(26) 



To investigate this point requires to find the expression 
of „„^ , and this is obtained by finding the expression 
of dp as a differential form with respect to d/S, d/3bar- In 
analogy with what previously done for other differential 
forms, we just define dp as the first order expansion of 
/AP) with respect to d/3, d/3bar- This gives, using 
again the Lemma, 

dp = -^{APAiH + pV)) d/3 + |(APA1/) d/3bar (27) 



from which one gets 



dp 



dP\ 



bar 



Uapav) . 



(28) 



On the other hand, if a suitable decorrelation occurs, one 
has as usual 



1 



{APAV)^{{P-{P)){V-{V))) 



(29) 



Finally, the r.h.s. of (|29p is just equal to — 1//3, as one 
finds in the Landau-Lifshitz textbook 11] . 

Thus after a sufficiently long time the temperature of 
the thermostat is an integrating denominator of the form 
5Q. It is presumable however that such a relaxation time 
might in general be rather short, especially in comparison 
with the relaxation times related to the specific heats at 
low temperatures. 



V. CONCLUSIONS 

In this paper we have shown how it is possible to for- 
mulate statistical thermodynamics for a Hamiltonian sys- 
tem, in terms of the processes induced by a change of the 
thermodynamic parameters. A characteristic feature of 
our approach is that it allows one to deal with metaequi- 
librium (or partial equilibrium) situations in which the 
time-correlations did not decay to zero, in addition to 
the equilibrium situations described in the famililar way 
through equilibrium ensembles. An approach somehow 
similar to ours was recently proposed by a C. Bernardin 
and S. OUa [iO|- The main difference is that we frame 
the problem in a time-independent Hamiltonian setting 
in which the processes are induced by nonequilibrium 
measures for the initial data, whereas in the quoted pa- 
per the processes are induced by imposing suitable time 
dependencies on the parameters of the global system. 

A relevant feature of the present approach is that it 
involves in a substantial way, in addition to the system 
of interest, also the measurement apparatuses. Conse- 
quently, one has the problem of understanding whether 
a physical, actually measured, property belongs to the 
system of interest or rather depends on the measuring 
procedure too. Standard equilibrium statistical thermo- 
dynamics, being formulated in terms of the Hamiltonian 
of the system of interest only, avoids in principle any ref- 
erence to the method of measurement, and perhaps one 



would be tempted to say that this is indeed correct. In 
general, however, things are not so. And actually, in the 
work [2, of which it was said that it gave rise to a new 
epoch of thermodynamics, one even finds the sentence: 
"How does one interpret a quantity such as the specific 
heat in a nonequilibrium situation? Clearly, the dynam- 
ics of how one performs the measurement enters into 
the measured value". This circumstance finds a counter- 
part in the present approach. Indeed the thermodynamic 
quantities we are considering, typically the specific heat 
of the system of interest, are in general represented by 
expressions which, through the intermediate of the mea- 
sure, involve not only the system of interest, but also the 
sources. So the measured quantities in general depend 
on the measurement procedure, and the dependence dis- 
appears only when a suitable relaxation has occurred, if 
indeed it takes place. 

However, the main interest of the proposal made here 
is that it may provide an answer to a question of principle 
that remains open within the equilibrium approach. The 
problem is to understand how many degrees of freedom 
should be attributed to a microscopic model. Why can 
one neglect the rotational degrees of freedom in a system 
of hard, completely smooth, spheres? Conversely, why 
does a supercooled liquid apparently have a number of 
degrees of freedom larger than the corresponding crystal? 
Problems of this type are usually met in chemistry or 
in glasses (see for example uM . In fact, however, they 
touch on questions of principle that affiicted very much 
Boltzmann (see |12, and the quotations reported in [l3|) 
and Gibbs [1J| , and were later taken up by authors such 
as Rayleigh and Jeans (see the first Solvay Conference 
15 1 , particularly the letter sent by Rayleigh) and Nernst 
16|, in connection with foundational problems of quan- 
tum mechanics. 

In an appendix it is discussed how the present results 
may be of interest for the Fermi-Pasta-Ulam problem, 
which also presents metaequilibrium situations (as ap- 
parently first suggested by Parisi and collaborators^), 
and thus constitutes a modern arena for dealing with 
foundational problems of the type just mentioned. 



Appendix 

Possible Relevance of Metaequilibrium Ther- 
modynamics for the Fermi Pasta Ulam Problem. 

In this paper we have shown how a statistical thermody- 
namics different from the equilibrium one can be formu- 
lated in situations in which the system does not present 
a full dynamical relaxation, i.e., in situations of partial 
equilibrium, as Landau and Lifshitz would say. This can 
occur even if the system was originally in an equilibrium 
state, described by an invariant Gibbs measure. In par- 
ticular, if the dynamics is such that the decorrelation of 
H + pV is negligible for long times when T — > 0, i.e., if 
one has ^H+pviT) — cr^.y for long times when T ^^ 0, 
then one has C„ ~ 0, and so the system behaves in agree- 



ment with the third principle. 

Situations of this type were conceived by Bohzmann 
since his celebrated paper published in the journal 
Nature |12| in the year 1895 (see also the quotations re- 
ported inll3h and were amply discussed in the year 1916 
by Walther Nernst|16|. who came to the point of con- 
ceiving that equipartition of energy might be compatible 
with Planck's law. According to Nernst, the key point 
consists in distinguishing between the mechanical energy 
(which corresponds to equipartition if the initial data are 
Gibbs distributed), and the energy that can be dynam- 
ically exchanged in a process. Only the exchangeable 
energy should be identified with the thermodynamic in- 
ternal energy U, while the complementary part should 
play the role of a kind of "zero-point energy". 

Foundational questions of this type were reconsidered 
later in the context of the FPU problem. As previously 
mentioned, the idea that in the FPU problem one meets 
with metaequilibrium situations analogous to those of 
glasses was first advanced by a group of people about 
Parisi in the year 1982 [4| (see also|^ and the review 0). 

In the FPU problem one has a classical approach, in 
which initial data are considered with only some of the 
normal modes excited (typically, the low frequency ones), 
and the result found is that, for low enough specific ener- 
gies, up to extremely long times energy remains confined 
among a certain group of low frequency modes. Already 
in the year 1972, in a paper of Cercignani et al. [17]| by 
the title "Zero-point energy in classical non-linear me- 
chanics" it was however proposed that even in the FPU 
problem one might meet with situations of the type con- 
ceived by Nernst. 

Numerical studies on the variant of the FPU problem 
with Gibbs distributed initial data came later, with re- 
sults that were variously interpreted (seell8l.[l9l). Finally 
a case was considered, which involves a situation very 
similar to the one discussed here (20|, where one mea- 
sures the specific heat of an FPU system (with an initial 
Gibbs distribution), put in contact with a heat reservoir. 
The result found is that at high enough energies a com- 
plete decorrelation occurs, and the standard equilibrium 
value is measured. Instead, at low enough energies situa- 
tions of an apparent metaequilibrium (or partial equilib- 
rium) are met, and smaller values of the specific heat are 
measured. Situations of metastability with hysteresis ef- 
fects had been previously reported [21| . Finally, the time 
autocorrelation function of the mode energies, or of the 
energies of packets of nearby frequencies, were numeri- 
cally studied in the paper [22|. 

So much for what concerns numerical results. In the 
meantime a consistent progress was achieved at an ana- 
lytical level, because it was proven that a weak version of 
perturbation theory can be formulated in the thermody- 
namic limit [23| . At variance with the classical approach, 
one renounces to control all initial data, and only looks 
at quantities such as the time-autocorrelation functions 
of suitable dynamical variables, estimated with respect 
to the Gibbs measure. This allows one to obtain results 



which hold in the thermodynamic limit of infinitely many 
particles, with a finite temperature. It was found that 
such autocorrelation functions do not decay to zero up to 
rather long times: a Nekhoroshev-like exponential stabil- 
ity was found for a <I>^-type model [2J|, while one step 
of perturbation theory was performed in the very FPU 
model '25i|. Such results appear to disprove the common 
expectation (see for example 10) that situations of par- 
tial equilibrium as those discussed here should not occur 
for infinite systems. 

Concerning the connections with the foundations of 
quantum mechanics, we now add a comment on the work 
of Poincare of the year 1912 on the necessity of quantiza- 
tion [2y| (see also the work[27l of Ehrenfest) . The problem 
was posed in the following way. Assume that the inter- 
nal energy U of the system of interest (the blackbody) be 
expressed as the mean value of the Hamiltonian H with 
respect to a suitable measure. 



U 



H{z)a{z)dz . 



(30) 



Then, it is asked which qualitative properties should the 
measure a have in order that the function U given by 
pop agree with the phenomenological internal energy. 
The answer of Poincare was that, under quite general 
assumptions, quantization is necessary, i.e., the measure 
has to be singular. 

Our comment is that, if one takes the metaequilibrium 
perspective, the claim that U should be expressed as an 
equilibrium mean value of the energy H (rather than a 
nonequilibrium mean value of A'^H), is unfounded. The 
result of Poincare shows that, if one does so, i.e., if one 
insists in representing U as an equilibrium mean value of 
_ff , then everything goes as if energy were quantized. Or, 
in Einstein's words at his contribution to the first Solvay 
Conference (seel28l, section 3), "The statistical properties 
of the phenomenon are the same as if energy were trans- 
ferred only by integer quanta of size hv ". 

Curiously enough, an opinion similar to the latter one 
of Einstein was expressed by Poincare too, in a paper of 
a less technical character [23| written shortly after the 
previously mentioned one, where he says: "Will disconti- 
nuity reign over the physical universe and will its triumph 
he definitive? Or rather will it he recognized that such a 
discontinuity is only an appearance and that it dissimu- 
lates a series of continuous processes? The first person 
that saw a collision helieved to be observing a discontinu- 
ous phenomenon, although we know today that the person 
was actually seeing the effect of very rapid changes of ve- 
locity, yet continuous ones. " And the conclusion was: 
"To try to express today an opinion about these problems 
would mean to be wasting one 's ink. " 

Now, our admiration for Poincare is unlimited, but our 
personal feeling, or rather hope, is that perhaps on this 
point Einstein was seeing farther than him. It appears 
that the main intuition of Einstein about the role of fluc- 
tuations was correct, although "rather ironically" (to use 
the words oi 3) he couldn't have a full appreciation of 



their dynamical features, which was obtained with FD 
theory only half a century later. It is just the stabiliza- 
tion of the time-autocorrelation functions to well definite 
nonvanishing values that allows for a thermodynamics 
different from the equilibrium one to hold. And this, not 
only for chemical systems or for glassy-like ones, but per- 
haps also for systems of FPU type, which are of interest 



for the foundations of physics. By the way, if one looks 
at the objections raised by Poincare againsts the "partial 
or incomplete equilibrium" point of view of Rayleigh and 
Jeans, one realizes that what he could not conceive is just 
that such partial equilibrium situations may actually oc- 
cur in nature, and that they may be dealt with through 
the methods of phenolemological thermodynamics. 



[1] L.D. Landau, E.M. Lifshitz, Statistical Physics, Pergamon 
Press (London, 1958). 

[2] N.O. Birge, Pliys. Rev. B 34, 1631 (1986); N.O. Birgc, S.R. 
Nagel, Phys. Rev. Lett. 54, 2674 (1985). 

[3] J.K. Nielsen, J.C. Dyre, Phys. Rev. B 54, 15754 (1996). 

[4] E. Fucito, F. Marchesoni, E. Marinari, G. Parisi, L. Peliti, 
S. RuflFo, A. Vulpiani, J. Phys. 43, 707 (1982). 

[5] G. Parisi, Europhys. Lett. 40, 357 (1997). 

[6] B. Benettin, A. Carati, L. Galgani, A. Giorgilli, in The 
FPU problem: a status report, G. Gallavotti ed.. Lecture 
Notes in Physics Vol. 728, Springer Verlag (Berlin, 2007). 

[7] E. Fermi, Thermodynamics, Dover (New York, 1956). 

[8] A.J. Khinchin, Mathematical Foundations of Statistical 
Mechanics, Dover (New York, 1949). 

[9] N.I. Chernov, J.L. Lebowitz, Ya.G. Sinai, Uspckhi Mat. 
Nauk 57 (6), 3 (2002). 
[10] C. Bernardin, S. OUa, Non-equilibrium macroscopic dy- 
namics of chains of anharmonic oscillators, preprint. 
[11] M. Mezard, G. Parisi, Phys. Rev. Lett. 82, 747 (1999). 
[12| L. Boltzmann, Nature 51, 413 (1895). 
[13[ A. Carati, L. Galgani, B. Pozzi, The problem of the rates 
of thermalization and the relations between classical and 
quantum mechanics, in M. Fabrizio et al. eds.. Mathemat- 
ical models and methods for smart materials, Advances 
in Mathematics n. 62, World Scientific (Singapore, 2002). 
Also available at Luigi Galgani's homepage. 
[14[ J.W. Gibbs, Elementary Principles in Statistical Mechan- 
ics, Dover (New York, 1960). 



[15[ La theorie du rayonnement et les quanta, P. Langevin and 

M. de Broglie eds., Gauthier-Villard (Paris, 1912). 
[16[ W. Nernst, Verhandl. Deutschen Phys. Gesell. 18, 83 

(1916). 
[17[ C. Cercignani, L. Galgani, A. Scotti, Phys. Lett. 38A, 403 

(1972); L. Galgani, A. Scotti, Phys. Rev. Lett. 28, 1173 

(1972). 
[18[ R. Livi, M. Pettini, S. Ruffo, A. Vulpiani, J. Stat. Phys. 

48, 539 (1987). 
[19[ A. Pcrronace, A. Tenenbaum, Phys. Rev. E 57, 100 (1998); 

G. Marcelh, A. Tenenbaum, Phys. Rev E 68, 041112 

(2003). 
[20[ A. Carati, L. Galgani, Europhys. Lett. 74, 528 (2006). 
[21| A. Carati, L. Galgani, J. Stat. Phys. 94, 859 (1999). 
[22| A. Carati, L. Galgani, A. GiorgiUi, S. Paleari, Phys. Rev. 

E 76, 022104 (2007). 
[23[ A. Carati. J. Stat. Phys. 128, 1057 (2007). 
[24| A. Carati, A. Maiocchi, Comni. Math. Phys. 314, 129 

(2012). 
[25[ D. Bambusi, A. Carati, A. Maiocchi, in preparation. 
[26| H. Poincare, J. Phys. Th. Appl. 2, 5 (1912); (Euvres IX, 

626. 
[27[ P. Ehrenfest Ann. der Phys. 36, 103 (1911). 
[28[ A. Einstein, L'etat actuel du probleme des chaleurs speci- 

fiques. Contribution to the first Solvay Conference. 
[29[ H. Poincare, Revue Scientifique 17, 225 (1912); (Euvres 

IX, 654. 



